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Abstract 


The Kudryashov method to look for the exact solutions of the nonlinear 
differential equations is presented. The Kudryashov method is applied to 
search for the exact solutions of the Liouville equation and the Sinh-Poisson 
equation. The equations of magnetohydrostatic equilibria for a plasma in a 
gravitational field are investigated analytically. An investigation of a family 
of isothermal magnetostatic atmospheres with one ignorable coordinate cor- 
responding to a uniform gravitational field in a plane geometry is carried out. 
The distributed current in the model J is directed along the x-axis where x 
is the horizontal ignorable coordinate. These equations transform to a single 
nonlinear elliptic equation for the magnetic vector potential u. This equation 
depends on an arbitrary function of u that must be specified. 


Keywords: Kudryashov method; magnetostatic equilibria; nonlinear evolu- 
tion equations; traveling waves. 


1 Introduction 


The equations of magnetostatic equilibria have been used extensively to 
model the solar magnetic structure [1,4,9,11]. An investigation of a fam- 
ily of isothermal magnetostatic atmospheres with one ignorable coordinate 
corresponding to a uniform gravitational field in a plane geometry is carried 
out. The force balance consists of the force between J A B (B, magnetic 
field induction, J is the electric current density), the gravitational force, and 
gas pressure gradient force. However, in many models, the temperature dis- 
tribution is specified a priori and direct reference to the energy equations is 


*Corresponding author 

Received 30 March 2015; revised 17 June 2015; accepted 1 August 2015 

N. Kadkhoda 

Department of Mathematics, Faculty of Basic Sciences, Bozorgmehr University Of Qaenat, 
Qaenat, Iran. e-mail: kadkhoda@buqaen.ac.ir 


H. Jafari 
Department of Mathematics and Computer Science,University of Mazandaran, Babolsar, 
Iran. e-mail: jafariQumz.ac.ir 


43 


44 N. Kadkhoda and H. Jafari 


eliminated. In solar physics, the equations of magnetostatic have been used 
to model diverse phenomena, such as the slow evolution stage of solar flares, 
or the magnetostatic support of prominences [20]. The nonlinear equilib- 
rium problem has been solved in several cases [3, 8,17, 18]. In this paper, 
we obtain the exact analytical solutions for the Liouville and sinh-Poisson 
equations using the Kudryashov method. Because these two models will be 
special cases of magnetostatic atmospheres model. Also here there is force 
balance between different forces. The Kudryashov method was developed by 
Kudryashov on the basis of a procedure analogous to the first step of the test 
for the Painlev property [2,7,7,9,10]. The paper is organized as follows : 

In Section 2, we describe the methodology of Kudryashov method for solv- 
ing nonlinear evolution equations when the Riccati equation is used as the 
simplest equation. We describe the Basic equations in Section 3. We apply 
this methodology and obtain exact solutions of the Liouville and sinh-Poisson 
equations in Section 4. Finally, the concluding remarks are presented in Sec- 
tion 5. 


2 Analysis of the Kudryashov method 


We consider a partial differential equation and we assume that by means of an 
appropriate transformation this partial differential equation is transformed 
to a nonlinear ordinary differential equation in the form 


Pa ot 0" a 0: (1) 


Exact solution of this equation can be constructed as finite series 
u(é) => Ala)’, (2) 


where G(€) is a solution of some ordinary differential equation referred to as 
the simplest equation. The simplest equation has two properties: 


1. the order of simplest equation should be less than the order of equation 
(1). 


2. we should know the general solution of the simplest equation or at least 
exact analytical particular solution(s) of the simplest equation. 


In this paper, we use the equation of Riccati, as the simplest equation. This 
equation is a well-known nonlinear ordinary differential equation which pos- 
sesses the exact solution constructed by elementary function. In this paper 
for the Riccati equation 


G'(€) = eG(é) + dG(g)’, (3) 
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we use the solution 


BOREL. 
OS) = TT aexplel + 01° 


d<0,c>0, (4) 


and 
c exp[c(€ + £0) ] 


1+ dexple(E+ £)]' 


Here € is a constant of integration. Now u(€) can be determined explicitly 
by using the following three steps: 


G(g) = 


d>0,c<0. (5) 


e Step (1). By considering the homogeneous balance between the highest 
nonlinear terms and the highest order derivatives of u(€) in equation 
(1), the positive integer n in (2) is determined. 


e Step (2). By substituting equation (2) with equation (3) into equation 
(1) and collecting all terms with the same powers of G together, the left 
hand side of equation (1) is converted into a polynomial. After setting 
each coefficient of this polynomial to zero, we obtain a set of algebraic 
equations in terms of A; (i = 0,1,2,...,n), c, d. 


e Step (3). Solving the system of algebraic equations and then substi- 
tuting the results and the general solutions of (4) or (5) into (2) gives 
solutions of (1). 


3 Basic equations 


The relevant of magnetohydrostatic equations consisting of the equilibrium 
equation with force balance will be as: 


JA B-pV®-VP=0, (6) 
which is coupled with Maxwells equations: 


VAB 
J= (7) 
Lb 


V- B=0, (8) 


where P, p ,jand ® are the gas pressure, the mass density, the magnetic 

permeability and the gravitational potential, respectively. It is assumed that 

the temperature is uniform in space and that the plasma is an ideal gas with 

equation of state p = pRoTo, where Ro is the gas constant and 7p is the 

temperature. Then the magnetic field B can be written by the following: 
Ou —Ou 


B=VuaA x B, ee (Beye : 
ul eg + Bye ( De By (9) 


46 N. Kadkhoda and H. Jafari 


The form of (9) for B ensures that V- B = 0, and there is no mono pole 
or defect structure. equation (6) requires the pressure and density be of the 
form [11]: 
= 1 =z 
Ply2)=Plwle®, ple) = rE PWwe®, (10) 


where h = Foto is the scale height. Substituting equations (7-10) in equation 


(6), we obtain 


V2utf(uje*® =0, 11) 
where ap 
fw) =e. 12) 
Equation (12) gives 
1 
P(u) = Po+ : / f(ujdu 13) 
Substituting equation (13) into equation (10), we obtain 
1 =z 
Ply.2)=(Po+= ff flujau)e*, (14) 
1 1 =z 
ply2)= (Pot = f Flu)du)e®, (15) 
g LM 


where Po is constant. Taking transformation 


a1 +itg=etet (16) 
equation (12) reduces to 
2 2 
care a +0? f(u) en? = 0, (17) 


Ont 048 


These equations have been given in Khater et al. (2000). 


4 Application of the Kudryashov method 


In this section, we will investigate the Kudryashov method for solving specific 
forms of f(u). 
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4.1 Liouville equation 


We first consider Liouville equation and the following equation will be special 
case of equation (17). Let us assume f(u) has the form (Dungey, 1953; Low, 
1975): 


f(u) = 0? Ape” 4, (18) 
where Ag and a? are constants. Hence 
2 A2 _2A 


P(y,z) = (Po + a AO ew, (19) 


Inserting equation (18) into equation (17) we obtain 


V2A/Ap = Pare ao ten), (20) 
where V? = 2 + a Let us set 
1 2 
A z 
As a + wy, 2), (21) 


where L is a constant. Then equation (20) becomes 


V2w — Pare 2b tH 72, (22) 
Let us identify 1 by 
2. 322 2 al 
ge One 2 
Re (23) 


and inserting equation (23) into equation (22) we obtain a Liouville type 
boa + ott coo “pe 7? =0. (24) 


In order to apply the Kudryashov method, we use the wave transformation 
€ =x — kt and change equation (24) into the form 


(1+ k*)¢" =a7le-*?, (25) 
we next use the transformation 
v=e*?, (26) 


we obtain a 
(1+ k?)uv" — (1 +k?) (v')? + 207 Pu? = 0, (27) 


with balancing according step (1) we get n = 2, therefore the solution of (27) 
can be expressed as follow: 
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v(§) =) A(G(e))’. (28) 


Substituting equation (28) along with (3) into (27) and setting the coefficients 
of all powers of G to zero, we obtain a system of nonlinear algebraic equations 
for Ag, Ay, Ag . Solving the resulting system with the help of mathematica, 
we have the following sets of solutions: 


Ao = 0, 

2 
A, =—- Sa, (29) 
Ag > a Crk ), 


where € = x—kt, A,a,l are constants. Therefore, substituting (29) in (28) and 
general solution (3) according to (4), we obtain solution of (27) as follows: 


c?d(1 +k?) e(e(E+€o)) 


m(6) = a2 (1 — de((E+€o)))2’ 30) 
whered<0,c>0,€=a-—kt. Using transformation 
v=e*?, (31) 
we get solution of (24) as follows: 
1 2d(1 2 (c(€+£o)) 
inf—£ (1+ k*) € (32) 


o1(€) = 9 [ PQ? (1 _ deere!” 


Now substituting (29) in (28) and general solution (3) according to (5), we 
obtain solution of (27) as follows: 


ed(1+k?) — elelé+éo)) 


v2(€) = Pa? (1 + dewréo)))2 (33) 
where d>0,c<0,€=2-—Akt. Using transformation 
poer?, (34) 
we get solution of (24) as follows: 
$e = ; we Es k2) 4 aaeangh (35) 
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4.2 Sinh-Poisson equation 


In this section, we consider sinh-Poisson equation which plays an important 
role in the soliton model with BPS bound [4,6]. Also, this equation will be 
special case of equation (17). If we assume 


2A 
flu) = — (2) sinh(9), (36) 
The same as above we have 
bua + Or = 8° sinh(¢). (37) 


In order to apply the Kudryashov method, we use the wave transformation 
€ =x — kt and change equation (37) into the form 


(1+ k2)4” = B? sinh(4), (38) 


we next use the transformation 


Vv = e?, 
ae = wes, oy 
we obtain 
2(1+ k?)ov" — 2(1 +k?) (v')? — 6?(v3 — v) =0. (40) 


With balancing according to step (1) we get n = 2, therefore the solution of 
(40) can be expressed as follows: 


v8) =D A(G(E))’. (41) 


Substituting equation (41) along with (3) into (40) and setting the coefficients 
of all powers of G to zero, we obtain a system of nonlinear algebraic equations 
for Ag, Ay, Ag . Solving the resulting system with the help of mathematica, 
we have the following sets of solutions: 


Ao = 1, 

Ay = 44 

Ay = 4 ou 
c= Tite ; 


where € = x — kt, 4,8 are constants. Therefore, using Substituting (42) in 
(41) and general solution (3) according to (4), we obtain solution of (40) as 
follows: 
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1 + dec(§+80))* 
(oo is d<0, c>0, (43) 
(1 = dec(§+80)) 
where c = oe When 6 > 0ore=- FA, when 6 < 0 and € = a2 — kt. 
Using transformation 
v=e?, (44) 


we get solution of (37) as follows: 


c(E-+€0))? 
1 ie + de ) 


(6) = lol ceaeny? (45) 


where c = Fer when 6 > O0Oorc = =a when 8 < 0. Now with 
Substituting (42) in (41) and general solution (3) according to (5), we obtain 


solution of (40) as follows: 


1 — decé+é)? 
vo(€ aie, d>0, c<0, (46) 
(1 + dec(E+£0)) 
where c = oe When 6 <0ore=- FA, when 8 > 0 and € = a — kt. 
Using transformation 
v=e?, (47) 


we get solution of (37) as follows: 


(ie dee(E+60))? 


= In] -————____ , ], 48 
aay (48) 
where c= 7555 when 8 <0 ore=— 755 when 6 > 0. 


5 Conclusion 


This study shows that the Kudryashov method is quite efficient and practical 
and is well suited for use in finding exact solutions for the Liouville and Sinh- 
Poisson equations. The reliability of the method and the reduction in the 
size of computational domain give this method a wider applicability. In this 
paper, the Kudryashov method has been successfully used to obtain some 
exact travelling wave solutions for the Liouville and Sinh-Poisson equations. 


Kudryashov method for exact solutions of isothermal ... 51 


Acknowledgements 


The authors are very grateful to the referees for their comments. 


References 


1. Heyvaerts, J., Larsy, J.M., Schatzman, M. and Witomsky, P. Blowing up 
of twodimensional magnetohydrostatic equilibria by an increase of electric 
current or pressure, Astronomy & Astrophysics 111 (1982) 104-112. 


2. Jafari, H., Kadkhoda, N. and Khalique, C.M. Travelling wave solutions of 
nonlinear evolution equations using the simplest equation method, Comp. 


Math. Appl 64 (2012) 2084-2088. 


3. Khater, A.H. Analytical solutions for some nonlinear two-dimensional 
magnetostatic equilibria, Astr. Space. Sci 162 (1989) 151-157. 


4. Khater, A.H., El-Attary, M.A., El-Sabbagh, M.F. and Callebaut D.K. 
Two-dimensional magnetohydrodynamic equilibria, Astr. Space. Sci 149 
(1988), 217-223. 


5. Khater, A.H., El-Kalaawy, O.H and Callebaut D.K. Backlund transfor- 
mations and exact solutions for a nonlinear elliptic equation modelling 
isothermal magentostatic atmosphere, IMA J. Appl. Math 65 (2000) 97- 
108. 


6. Kudryashov, N.A. One method for finding exact solutions of nonlinear dif- 
ferential equations, Commun. Nonlin. Sci. Numer. SimuL 17 (2012) 2248- 
2253. 


7. Kudryashov, N.A. Simplest equation method to look for exact solutions of 
nonlinear differential equations, Chaos. Solit. Fract 24 (2005) 1217-1231. 


8. Lerche, I. and Low, B.C. On the equilibrium of a cylindrical plasma sup- 
ported horizontally by magnetic elds in uniform gravity, Solar. Phys 67 
(1980) 229-243. 


9. Lerche, I. and Low, B.C. Some nonlinear problems in astrophysics, Physica 
D: Nonlinear Phenomena 4 (1982) 293-318. 


10. Low, B.C. Nonlinear Force-Free Magnetic Fields, Rev. Geophys. Sp. Phys 
20 (1982) 145-159. 


11. Low, B.C. Evolving force-free magnetic fields. III- Stages of nonequilib- 
rium and the preflare stage, The Astrophys. J 239 (1980) 377-388. 


D2 N. Kadkhoda and H. Jafari 


12. Martinov, N.k. and Vitanov, N.k. On the correspondence between the self- 
consistent 2D Poisson- Boltzmann structures and the sine-Gordon waves, 


J. Phys. A: Math. Gen. 25 (1992) L51. 


13. Martinov, N.K. and Vitanov, N.K. On self-consistent thermal equilib- 
rium structures in two-dimensional negative-temperature systems, Cana- 


dian Journal of Physics 72 (1994) 618-624. 


14. Salehpour, E., Jafari, H. and Kadkhoda, N. Application of (Se 
expansion method to nonlinear Lienard equation, Indian Journal of Sci. 
and Tech 5(2012) 2454-2456 . 


15. Vitanov, N.K., Dimitrova, Z.I. and Kantz, H. Modified method of simplest 
equation and its application to nonlinear PDEs, Appl. Math. Comput. 216 
(2010) 2587-2595. 


16. Vitanov, N.K. Modified method of simplest equation: Powerful tool for ob- 
taining exact and approximate traveling-wave solutions of nonlinear PDEs, 
Commun. Nonlin. Sci. Numer. SimuL 16 (2011) 1176-1185. 


17. Webb, G.M. Isothermal magnetostatic atmospheres. II - Similarity solu- 
tions with current proportional to the magnetic potential cubed, Astrophys. 
J 327 (1988) 933-949. 


18. Webb, G.M. and Zank, G.P. Application of the sine-Poisson equation in 
solar magnetostatics, Solar. Phys 127 (1990) 229-252. 


19. Wazwaz, A.M. Partial Differential Equations and Solitary Waves Theory, 
Springer Verlag, (2009). 


20. Zwingmann, W. Theoretical study of onset conditions for solar eruptive 
processes, Solar. Phys 111 (1987) 309-331. 


Led (ob gevcolinn pinweil G59 cle Gl ym gly doles oy 5 role sy 


Y epee peer \ is aS all Cae 


els 298 al pole Stile LIB gS) 59 olStsls : 
Pagel gle 9 ely osKtsls cghsjb lRssls r 


VAY slay Ve lie dak VAP sb, 5 Vaud ciel dlis cibys VTA ial 18 dle Gib,» 


veel od SI! ed ne ntl io Vales gs cl lye opener oly USL oS da) 2 ode 
DS Opell role grew doles 9 vagd doles Gibs clelse grrr cly GdL oS Ss) 
Eig da BETIS lnc Ss 39 Vahl del yas ollie Realy ly ySta IslatiesMlaley Waals ty 


Slate Sb bes pe pebbles nail bol gh S guy wl ob 15 Quy ory chbs 
9 J Jue 59 fewer ws Oby> Sane ri ter ples! dbs durio 50 SS1gS HLS line SG 4 ble 
doles So 4 und Vales cyl cael ond Cuban ahh (ce Hd) Slits X Ol 59 SX pyres slucel 
Kosly U alyrds OE Sy 4 doles oul gd ge U ecbbite Gila y Gentle oly ed nb pay 

195 adeeb S Cul 


KS 5 Clael bd st AIS Vo lee tb lite Jal abby S say + GAT OLUS 


